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' Abstract: In this paper, we examine the interacting dark energy model in f(T) cosmology. 

We assume dark energy as a perfect fluid and choose a specific cosmologically viable form 

I ^ . f{T) = PVT. We show that there is one attractor solution to the dynamical equation of 

f(T) Friedmann equations. Further we investigate the stability in phase space for a general 

f(T) model with two interacting fluids. By studying the local stability near the critical 

Q ' points, we show that the critical points lie on the sheet u* = (c — l)v* in the phase space, 

I ■ spanned by coordinates (u, v,Q, T). From this critical sheet, we conclude that the coupling 

5h . 

between the dark energy and matter c € (—2, 0). 
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I. INTRODUCTION 

(N ■ 
t-H ■ 

Astrophysical observations indicate that nearly seventy percent of the cosmic energy density is 
hidden in some unknown 'dark' sector commonly called as 'dark energy' (DE) Q,Q{. The remaining 
contribution to the total energy density is contained in dark matter and meager baryonic matter 
[jj. Dark energy is described phenomenologically by an equation of state (EoS) pd = wpd, where 
Pd and pd are the pressure and energy density of dark energy, while w is an EoS parameter which 
effectively describes cosmic acceleration. To describe DE, one must have pd < and w < 0. 
General relativity (GR) offers only one solution to this puzzle, , namely 'cosmological 
constant' which suffers from fine tuning and coincidence problems. Since GR fails to explain the 
cosmic accelerated expansion, one needs to modify curvature or the matter part in the Einstein 
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[4], scalar-tensor gravity 



field equations. Some notable examples are f(R) gravity [J], scalar-tensor gravity [5[ and Lovelock 
gravity f| and more recently f(T) gravity [71, to name a few. Other models where matter action 
is modified include the bulk viscous stress 8] and the anisotropic stress Q or some exotic fluid 
like Chaplygin gas (CG) [101 ] . One of the crucial tests to check the viability of extended 
theories of gravity is the potential detection of gravitational waves [11]. 

We assume a phenomenological form of interaction between matter and dark energy since these 
are the dominant components of the cosmic composition, following fl^ . The exact nature of this 
interaction is beyond the scope of the paper and unresolved till we obtain a consistent theory of 
quantum gravity. In these interacting dark energy-dark matter models, the dark energy decays 
into matter at a rate proportional to Hubble length. The interacting dark energy scenario 
can successfully resolve the coincidence problem then stable attractor solutions of the 
Friedmann- Robertson- Walker (FRW) equations can be obtained 1 131. Some observa- 

nn 

tional support to these models come from the astrophysical observations [14. 1151]. For 
motivation from field theory and p arti cle physics of the interacting dark energy, the 
interested reader is referred to [la . \v\ . 

The present paper is devoted to the study of dynamics of interacting dark energy in f(T) 
cosmology. This theory is based on torsion scalar T rather than curvature scalar R. We assume 
the dark energy in the form of a perfect fluid interacting with the matter. We study this model 
by the local stability method and than extend it for dark energy satisfying a more 
general form of equation of state. 



II. BASIC EQUATIONS 

A. Basics of f(T) gravity 

General relativity is a gauge theory of the gravitational field. It is based on the equivalence 
principle. However it is not necessary to work with Riemannian manifolds. There are 
some extended theories such as Riemann-Cartan, in them the geometrical structure 
of the theory has non-vanishing object of non-metricity. In these extensions, there are 
more than one dynamical quantity (metric) . For example this theory may be constructed from the 
metric, non-metricity and torsion [181 ]. Ignoring from the non-metricity of the theory, we 
can leave the Riemannian manifold and go to Weitznbock spacetime, with torsion and 
zero local Riemann tensor. One sample of such theories is called teleparallel gravity in which 
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we are working in a non-Riemannian manifold. The dynamics of the metric is determined 
using the scalar torsion T. The fundamental quantities in teleparallel theory are the 
vierbein (tetrad) basis vectors e* . This basis is an orthogonal, coordinate free basis, defined 
by the following equation 



This tetrad basis must be orthonormal and rjij is the Minkowski metric. It means that e^e w - = 5j. 
There is a simple extension of the teleparallel gravity, which is called f(T) gravity. 
In this theories, f is an arbitrary function of the torsion T. One suitable form of action for f(T) 
gravity in Weitzenbock spacetime is [3] 



S = J^I d4x€{T + /(T) + Lr ' 



Here e = dei(e* ), k 2 = 8irG. The dynamical quantity of the model is the scalar torsion T and L m 
is the matter Lagrangian. The field equation can be derived from the action by varying 



the action with respect to 



e-%(eS t H(l + h) ~ e t x T P uX S p ^f T 



iiX p 

1 

V 



+S, ^d,(T)f TT - -e\(l + f(T)) = ^GefT/, (1) 



where T v is the energy- momentum tensor for matter sector of the Lagrangian L m , defined by 

2 S I [V-9 L md 4 x] 



T, 



Here T is defined by 



where 



Sg 



rri Q f_LUrpp 



s p ^ = 1(k^ p + s^t%-s;t%), 



and the contorsion tensor reads K^ v p as 
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It is straightforward to show that this equation of motion reduces to Einstein gravity 
when f(T) = 0. Indeed, it is the equivalency between the teleparallel theory and 



Einstein gravity 20]. The theory has been found to address the issue of cosmic acceleration in 
the early and late evolution of universe [2l{] but this crucially depends on the choice of suitable 
.f(T), for instance exponential form containing T cannot lead to phantom crossing 



22]. Reconstruction of f(T) models has been reported in [231 ] while thermodynamics of f(T) 



cosmology including the generalized second law of thermodynamics has been recently 
investigated 



B. Interacting dark energy in f(T) cosmology 

We adopt the metric in the form of a flat Friedmann-Lemaitre-Robertson- Walker metric with 
metric ds 2 = dt 2 — a(t) 2 (dx l dxi), i = 1,2,3. We start with the Friedmann equation for the 
f(T) model Q 



1 + 2fr v 3 r 6 

where p = p m + Pd, and p m , pd represent the energy densities of matter and dark energy. 
The second FRW equation is 

-v. 

+ fr + 2T/tt 

For a spatially flat universe (k = 0), the total energy conservation equation is 



H ~ 2 (l+fr + 2Tf TT )- (3) 



p + 3H(p + p) = 0, (4) 

where H is the Hubble parameter, p is the total energy density and p is the total pressure of the 
background fluid. The so-called energy-balance equations corresponding to dark energy and dark 
matter are 25] 



p d + 3H(p d + p d ) = -Q, 

p m + 3Hp m = Q, (5) 

Here Q is the interaction term which corresponds to energy exchange between dark energy and 
dark matter. The function Q has dependencies on the energy densities of the dark matter and dark 
energy and the Hubble parameter, i.e. Q(Hp m ), Q(Hpd) or Q(Hp m , Hpd) 26]. Since the nature 



of both dark energy and dark matter is unknown, it is not possible to derive T from first principles. 
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To give a reasonable Q, we may expand like Q(Hp m , Hpd) ~ a m H p m + a^H p^. Since the coupling 
strength is also not known, we may adopt just one parameter for our convenience; hence we take 
ot m = OLd = c [273] . We here choose the following coupling function Q = 3Hc(pd + p m ). 

To perform the stability analysis of the cosmological model, it is always convenient to define 
dimensionless density parameters via 

K 2 Pd K 2 p d K 2 p m 

U= 3i^' US afl*' = MP' (6) 

Moreover the equation of state parameter of dark energy 

w = — = -. 7 

Pd u 

Also the equation of state of all the fluids together is 

- P d V (Q\ 

Wtot = ; = ~ pp. 8) 

Pd + Pm 1 - 

Using 0, we can rewrite ([2]) in dimensionless form 

n m = i-u-n T , (9) 

where Qt = ^p- — 2/t, is another dimensionless density parameter constructed for torsion scalar. 
In f(T) and teleparallel gravities, phase space analysis for different dark energy models 



has been reported in 



28, 



29] 



III. STABILITY ANALYSIS OF PERFECT FLUID 



In this section, we treat dark energy as a perfect fluid. Physically it means that the fluid has 
no viscosity or heat transfer property. Its a simple fluid which cannot be self-gravitating and has a 
smooth distribution in space. A perfect fluid is represented by a phenomenological linear equation 
of state connecting energy density and pressure by 

Pd = wp d , (10) 

where w is a 'constant' of proportionality but can be a function of time or e-folding parameter 
x = In a. For a general dark energy paradigm, the minimum condition to be satisfied is w < —1/3. 
This opens a window to construct variety of theoretical models to explain cosmic 
acceleration. 

Some of these well-known models are cosmological constant (w = —1), quintessence (w < —1/3) 
and phantom energy (w < —1). 
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The dynamical system representing the dynamics of perfect fluid's density and pressure reads 



as 



^ = -3[u + i; + c(l-n T )] 
ax 



l-tt T +v 



1 + 2Tf T + f TT 
j-u + cfl- 

i - n T + v 



+3u 



dv 

— = -3w[u + v + c(l -Qt)] 
ax 



-3v 



(11) 



(12) 



.l + 2Tf T + f TT . 

We notice that our dynamical system is under-determined: three unknown functions (u, v, Qt) and 



two coupled equations. Thus we must choose f(T) to solve the system. It is discussed in [30(] that a 
power-law form of correction term f(T) ~ T n (n > 1) such as T 2 can remove the finite-time future 
singularity. However, when n = 0, the correction term behaves like a cosmological constant. The 
model with n = 1/2 can be helpful in realizing power-law inflation, and also describes 
little-rip and pseudo-rip cosmology [30]. Due to these reasons, we choose 

f(T) = PV^T, 

where /3 is a constant. Note that choosing /3 = (or f(T) = 0) leads to teleparallel gravity 
which is equivalent to General Relativity. Note that this choice f(T) = has 
correspondence with the cosmological constant EoS in f(T) gravity [31] . This f(T) 
model can be recovered via reconstruction scheme of holographic dark energy [35]. 
Als o it can be inspired from a model for dark energy model form the Veneziano ghost 



32]. Recently Capozziello et al [33] investigated the cosmography of f(T) cosmology 
by using data of BAO, Supernovae la and WMAP. Following their interesting results, 
we notice that if we choose f3 = \f3/2Ho (ri m o — 1), than one can estimate the parame- 
ters of this f(T) model as a function of Hubble parameter Hq and the cosmographic 
parameters and the value of matter density parameter. 

Recently attractor solutions for the dynamical system with three fluids (dark matter, dark en- 
ergy and radiation) interacting non-gravitationally have been investigated to resolve the coincidence 
problem using similar f(T) 29]. 

To perform the stability of the system comprising equations (llip and (I12p . we first calculate 



the critical points (w*,v*) by equating 



ax ax 
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Point 


Wtot = V* 


w = — 


q = -1 - §(w* + «♦) 


A 


-2 


2 

c-1 


3c- 1 


B 


28 


W 


l-3*(l + £) 



TABLE I: Cosmological implications. 



We find two critical points A and B for system (jlip and (|12p given in Table-I. We check the stability 
of the dynamical system in the neighborhood of these critical points u = u* + Su, v = v* + Sv. This 
is performed by linearizing the system of equations of motion for u and v like 



d5u 3 . n ( -i 1 i 
— — = -u*ou + 31—1 + -u* ] dv, 



dSv 
dx 



-3w5u + 3 (1 — w + u*) 



(13) 
(14) 



Next we construct a Jacobian matrix from the coefficients of Su and Sv in the linearized (or 
perturbed) system and finding the eigenvalues from it. For the system (jlip and ()12p . the Jacobian 
matrix for any critical point is 



J 



3 ~\~ 2 



-3tu — 3tt; + 3 + 3^* 



(15) 



We find the eigenvalues Ai t 2 of the Jacobian matrix for the two critical points given in Table-I. 
A critical point (also called equilibrium point) is said to be stable if the corresponding eigenvalues 
are negative Ai^ < for all the values of the model parameters. Such stable critical points are 
called attractor solutions of the dynamical system i.e. solution of the dynamical system for 
different initial conditions which converge to a stable critical point. A critical point is said 
to be unstable if both the eigenvalues are positive, and a saddle point if one of the eigenvalues 
is positive. We are interested only in stable critical points. Sometimes we get conditionally stable 
critical points i.e. the critical point which can be stable only under some conditions on the model 
parameters. In cosmology conditionally stable points are also of interest. In Table-IH we write 
down the values of relevant cosmological parameters including the total equation of state and the 
deceleration parameter. Note that for f(T) = pVT, we have = 0, therefore wtot = v*- For A, 
q < when c < 1/3. 



We observe that A is an unstable critical point and B is conditionally stable. This means 
that if the interaction parameter c < this point is a stable critical point. Since our two 
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Point 




Ai 


A 2 


Stability Condition 


A 


(2(1 -c), -2) 


3(-l - § + ±Vw 2 +4wc) 


3(-l - f - ±Vw 2 +4wc) 


Unstable point 


B 




3(1 + 6) 


3(5 -w) 


Stable point if c < ^ 



TABLE II: Critical points and stability of the perfect fluid model ( 8 = ^(w — y/w 2 + Awe), w 2 +4wc> 
)■ 



dimensional phase space embedded in a three dimensional space, it has no chaotic behavior, 
and we can investigate the dynamical behavior of the the system by using the usual 
dynamical systems approach. In figure 1, we plotted the dynamical phase space of the (jlip . 
(|12p for different values of the parameters and for a range of the initial conditions. We adopted 
the initial value (here initial value means the present value) for energy density of dark energy 
u(0) = 0.7, but different values of pressure parameter -u(O) for different forms of dark energy 
in order to have a good agreement with the observational data. From figure-1 it turns 
out the trajectory of phantom energy in phase space is more steep than quintessence 
and cosmological constant. From figure-2, we plot the deceleration parameter for different dark 
energy models in f(T) cosmology. For cosmological constant, the trajectory starts near x = — 1 
(or a = 0.36) and terminates close to x = 4 (or a = 54.59), meanwhile the deceleration parameter 
remains fixed between x = —1 and x = 4 since Hubble parameter is constant. For quintessence, 
the deceleration parameter begins increasing from x = —0.5 (o = 0.6) to x = 1.5 (a = 4.48) where 
q vanishes. Thus deceleration q > is not possible with quintessence dark energy. For phantom 
energy, q starts decreasing from x = —1.3 but gets stable to value q = —0.7. 

IV. THE PHASE SPACE OF THE INTERACTING DARK ENERGY MODELS: 

GENERAL CASE 

For a generic form of the f{T), the system of equations (jlip . (|12p is not closed i.e. a third 
differential equation for is needed. The EoS of the fluid is in the form = $(/Orf). The general 
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-7- 

2 4 6 8 10 



FIG. 1: Phase space for perfect fluid form of dark energy. Three forms of dark energy are shown in figure: 
phantom energy v(0) = —0.84 (w = —1.2, green); cosmological constant v(0) = —0.7 (w = —1, blue); 
quintessence v(0) = —0.231 (w = —0.33, red). 
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FIG. 2: Behavior of deceleration parameter for perfect fluid form of dark energy. Three forms of dark energy 
are shown in figure: phantom energy v(Q) = —0.84 (w = —1.2, green); cosmological constant v(Q) = —0.7 
(w = —1, blue); quintessence v(0) = —0.231 (w = —0.33, red). 



dynamical system then reads 

du 
dx 



dv 
dx 



-3[c(l - fi T ) + u + v] 
3n(l - n T + v) 

+ 1 + / T + 2T/ TT ' liDj 
-3&[c(l-n T ) + u + v] 

+ Mi -Qt + v), (17) 



1 + /t + 2T/ tt ' 
cMt = 3T(1 - n T + v)(Tf T -f- 2T 2 f TT 
dx 1 + f T + 2Tf TT 



(18) 
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Here $' = The system (|16H18p is non- autonomous due to presence of / and we need to add 
another equation to it. From the definition of the T = —6H 2 it is easy to show that the fourth 
equation is 

3T(i - n T + v) 



dx 



(19) 



1 + fr + 2Tf TT 

Now the system (|16H19p is closed. We discussed the different possible attractors of the system 
$MM in Table-HE 



TABLE III: Critical Points of the general f(T) model. 



Critical sheet 


(u*,v*,Sl* T ,T*) 


Stability Condition 


C 
D 


(u*,v*,n* T (0),0) 
((c-l)v*,v*,T*,,n* T ) 


Physically unacceptable 
Conditionally stable 



Case C is not physically viable since when T* = it implies H = 0. It means that the local 
geometry in the neighbourhood of this point is Minkowski flat. We are not interested to such 
cases since astrophysically H ^ 0. But for case D, we conclude that the critical points in the phase 
space lie on the two dimensional surface 



u 



(c-l)v*. 



(20) 



Since EoS is w = this equation tells us that the critical point can be written in the form 

^Tj-. Since —1 < w* < — g, it shows that the i nter acting coupling 



■f^r. It shows that w* 



must be in the range c 6 (—2, 0), which is also reported in an earlier work 36]. This is a new 
constraint on the coupling constant c from the phase analysis approach. But from this analysis we 
can not obtain any new information about the values of the T*, Wp. We remember here that, this 
discussion is free from any specific form of the f(T). We conclude that for interacting dark energy 
models in f(T) gravity, there is only one physically acceptable critical point D. 



V. CONCLUSION 



In this paper, we investigated the stability and phase space description of a perfect 
fluid form of the dark energy interacting with matter. We wrote the general dynamical 
system equations for two fluids and found the critical points. We linearized the system of 
equations and showed that for perfect fluid case there is only one attractor solution. 

For any value of the EoS parameter w, we obtained a new bound on the coupling constant c < 



11 



It is a new constraint on c in the context of the f(T) gravity which shows that energy is being 
transferred from matter to dark energy. Thus we obtain a scenario where the universe becomes 
increasingly dark energy dominated. Further for a general f(T) model, there is only one critical 
point, and this critical point lives on the surface with equation u* = (c — l)v* in the four 
dimensional phase space which is spanned by four coordinates X = (u*, v*, Oy, T*). 
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